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Abstract 

The volume fluctuations in statistical mechanics are discussed. First, the volume fluctuations in 
ensembles with a fixed external pressure, the so called pressure ensembles, are considered. Second, 
a generalization of the pressure ensembles is suggested. Namely, the statistical ensembles with the 
volume fluctuating according to externally given distributions are considered. Several examples and 
possible applications in statistical models of hadron production are discussed. 
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I. INTRODUCTION 



Successful application of the statistical model to description of mean hadron multiplicities 
in high energy collisions (see, e.g., recent papers [l| and references therein) has stimulated 
investigations of properties of statistical ensembles. Whenever possible, one prefers to use the 
grand canonical ensemble (GCE) due to its mathematical convenience. The canonical ensemble 

h 

(CE) [2j should be applied when the number of carriers of a conserved charges is small (of 
the order of 1), such as strange hadrons sj], anti-baryons Q], or charmed hadrons Q]. The 
micro-canonical ensemble (MCE) 6| has been used to describe small systems with fixed energy, 
e.g. mean hadron multiplicities in proton-antiproton annihilation at rest. In all these cases, 
calculations performed in different statistical ensembles yield different results. This happens 
because the systems are 'small' and they are 'far away' from the thermodynamic limit (TL). 
The mean multiplicity of hadrons in relativistic heavy ion collisions ranges from 10 2 to 10 4 , 
and mean multiplicities (of light hadrons) obtained within GCE, CE, and MCE approach each 
other. One refers here to the thermodynamical equivalence of statistical ensembles and uses 
the GCE for calculating the hadron yields. 

Measurements of a hadron multiplicity distribution P(N) in interactions, including nucleus- 
nucleus collisions, open a new field of applications of the statistical models. The particle 
multiplicity fluctuations are usually quantified by the ratio of variance to mean value of a 
multiplicity distribution P(N), the scaled variance, and are a subject of current experimental 



activities. In statistical models there is a qualitative difference in the properties o 



mean 



multiplicity and scaled variance of multiplicity distributions. It was recently found 
11 , H . [13I that even in the TL corresponding results for the scaled variance are different in 



different ensembles. Hence the equivalence of ensembles holds for mean values in the TL, but 
does not extend to fluctuations. 

A statistical system is characterized by the extensive quantities: volume V, energy E, and 
conserved charge(s) 1 Q. The MCE is defined by the postulate that all micro-states with given 
V, E, and Q have equal probabilities of being realized. This is the basic postulate of the 
statistical mechanics. The MCE partition function just calculates the number of microscopic 



1 In statistical description of the hadron or quark-gluon systems, these conserved charges are usually the net 
baryon number, strangeness, and electric charge. In non-relativistic statistical mechanics, the number of 
particles plays the role of a conserved 'charge'. 
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states with given fixed (V, E, Q) values. In the CE the energy exchange between the considered 
system and 'infinite thermal bath' is assumed. Consequently, a new parameter, temperature 
T is introduced. To define the GCE, one makes a similar construction for conserved charge Q. 
An 'infinite chemical bath' and the chemical potential \x are introduced. The CE introduces 
the energy fluctuations. In the GCE, there are additionally the charge fluctuations. The 
MCE, CE, and GCE are the most familiar statistical ensembles. In several textbooks (see, 
e.g., Ref. Q,H|), the pressure (or isobaric) canonical ensemble has been also discussed. The 
'infinite bath of the fixed external pressure' po is then introduced. This leads to the volume 
fluctuations around the average value. 

A more general concept of the statistical ensembles was suggested in Ref. [16|. The statistical 
ensemble is defined by an externally given distribution of extensive quantities, P a (A). The 
construction of distribution of any property O in such an ensemble proceeds in two steps. 
Firstly, the MCE O-distribution, P mce (0; A), is calculated at fixed values of the extensive 
quantities A = (V, E, Q). Secondly, this result is averaged over the external distribution P a (A) 



P a (0) = J dAP a (A) P mce (0;A) . (1) 

The ensemble defined by Eq. ([I]), the a-ensemble, includes the standard statistical ensembles 
as particular cases. 

Recently, the micro-canonical ensemble with the volume fluctuations was introduced 13] for 
modelling the hadron production in high energy interactions. An introduction of the volume 
fluctuations was necessary in order to reproduce the KNO-scaling 0] of the hadron multiplicity 
distribution observed in proton-proton and proton-antiproton collisions at high energies. The 
volume fluctuations lead also to the power law shape of the single particle spectra at high 
transverse momenta. 

In the present paper the statistical ensembles with volume fluctuations are studied. 
First, the pressure ensembles are considered. In general, there are 3 pairs 2 of variables - 



2 In the present study we do not discuss the role of the total 3-momentum. As shown in Ref. [l3| the total 
momentum conservation is not important in the TL for thermodynamical functions and fluctuations in the 
full phase space. It may however influence the particle number fluctuations in the limited segments of the 
phase space. 
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{V,po), (E,T), (Q,fi) - and, thus, the 8 statistical ensembles 3 can be constructed. Among 
these 8 ensembles there are 4 pressure ensembles: (po,E,Q), (p ,T,Q), (p ,E,fi), and 
(po,T,/j,). In addition to the pressure canonical ensemble known from the literature, three 
other possibilities - pressure micro-canonical, pressure grand micro-canonical, and pressure 
grand canonical ensembles - are constructed and studied. In Section II, non-relativistic sta- 
tistical systems are discussed, whereas Section III presents the results for an ultra-relativistic 
ideal gas. 

In Section IV the concept of pressure ensembles is extended to the case of more general 
volume fluctuations. Namely, the statistical ensembles with volume fluctuating according to 
externally given distributions are introduced. The Summary presented in Section V closes the 
paper. 

II. NON-RELATIVISTIC BOLTZMANN GAS 
A. Canonical and Micro- Canonical Ensembles 

In this Section the system of non-relativistic Boltzmann particles is discussed. The (V, T, N) 



Canonical Ensemble (CE) partition function of N particles reads 19]: 



ZUV,T,N) = - J ... exp 



S(xi,...,xjv;Pi,...,Pjv) 



(2) 



T 

where T is the system temperature, the particle degeneracy factor is assumed to be equal to 
1, and E is the microscopic iV-particle energy usually presented as the sum of potential and 
kinetic terms, 

N p 2 

S(x 1 ,...,x iV ;pi,...,pjv) = C/(xi,...,xjv) +$^2^' ( 3 ) 

i=l 

with m being the particle mass. For the iV-particle energy given by ([3]) the integration over 
momentum in Eq. (j2j) can be done explicitly, 



(4) 



3 For several conserved charges {Qi} the number of possible ensembles is larger, as each charge can be treated 
either canonically or grand canonically. 
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The particle coordinate xi, . . . , are integrated over the system volume V. The CE thermo- 
dynamical functions can be expressed in terms of the free energy, 



F(V,T,N) = — T In Z ce {V,T,N) . 



(5) 



The CE pressure, entropy, chemical potential, and average energy are: 

dN 



V 



S 



-(-) 



E 



F + TS . (6) 



T,V 



For the non-interacting particles, the potential energy vanishes, f/(x 1; . . . , X;y) = 0, thus, the 
ideal gas free energy is: 



F id (V,T,N) 



NT - NT In 



V fmT\ 3/2 ' 
N \2n) 



(7) 



where we have assumed N ^> 1 and, thus, lniV! = N In N — N. The thermodynamical functions 
of the ideal gas from Eqs. (!6lT7|) read: 



P=^, S = ^N + N\n 



V_ /mT\ 3/2 
N \~2n) 



-Tin 



V /mT\ 3/2 ' 
N 1 2tt J 



E = -NT . 
2 



The (V, E, N) Micro-Canonical Ensemble (MCE) partition function of N non-relativistic 
Boltzmann particles reads: 

N 9 



If dx 1 dp 1 dx N dp 
Z mce (V,E,N) = - J —— ... — 

The MCE entropy is defined as, 



A' 
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E- C/(xi,...,xjv) 



Z—' 2m 

3=1 



(9) 



S(V,E,N) = ln[E Z mce {V,E,N)\ 



(10) 



where Eq is an arbitrary constant with a dimension of energy. For non-interacting particles 
U (xi, . . . , xat) = 0, one finds, 

V ( mE\ 3/2 ~ 
N \3Nn) 



V N (Em I2tt^ 3N I 2 ^ 
Z mce (V, E,N)= * { r ' I , r , S(V, E, N) = -N + Nln 



(11) 



It has been assumed that N > 1 and E > E , thus, r(3iV/2) = 3N/2 [ln(3iV/2) - 1], 
N\ = N(\n N — 1), and \n(E ■ E 3N ' 2 - 1 ) = \n(E 3N / 2 ). 

The MCE temperature, pressure, and chemical potential are the following: 



For non-interacting particles they read, 



1 3 N p N [i 

T ~ 2~E ' T ~ V 1 T 

If E — E, the ideal gas system of N particles in the volume V has the same temperature, 
pressure, and entropy in the MCE (I11H13P and in the CE ([8]). This means the thermodynamical 
equivalence of CE and MCE at N ^> 1. 



N \3Ntt ) 



3/2 



(13) 



B. Pressure Canonical and Pressure Micro-Canonical Ensembles 

The particle coordinates and momenta for the considered system are denoted as xi, . . . , Xjv 
and pi, . . . , p^, respectively. The E given by Eq. d3J) denotes the system energy, and V is the 
system volume. The 'thermostat' is now introduced with corresponding particle coordinates, 
Xi, Xat t , and momenta, Pi, Pat t , and with Et and Vr being the thermostat 

energy and volume, respectively. The system plus thermostat is described by the MCE, i.e. 
the total energy and the total volume are assumed to be fixed: E + Et = E* = const , 
V +Vt = V* = const . The probability distribution of system particle coordinates and momenta 
is proportional to, 



„ n t / N T p 2 

Mx 1; x^; Pl , PN ) ex HdPXidPPiSlE + Y.^k 

J i=i V j=i 



E* ) , (14) 



where the thermostat particle coordinates Xi,...,Xjv t are integrated over the thermostat 
volume Vt, and the system particle coordinates xi, . . . , x^ are contained in the volume V. The 
particles in the thermostat are assumed to be non- interacting. Thus, E T = Eji P i/( 2M ), 
with M being the particle mass. The momentum integration in Eq. (|T4"|) gives, 

N T p2 \ / rp \ 3N T /2 -1 



1 



( T P 2 \ / r? \ ' 

E+ Hw - E j K ^ - e)3Nt/2 ~ 1 ^i 1 '^) 



E \ 3W T /2 -1 , E N 3Af T /2 



1 - E~ T ) S I 1 - 3i^372.» " 6XP( - B/T) • (15) 



where it has been assumed, N T — > oo, E/E T — > 0, and the ideal gas equation for the thermostat 
energy, E T = 3TN T /2, has been used. The integration over thermostat particle coordinates X, ; 
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in Eq. (JHj) gives the factor depending on the system volume V, 



d 3 Xi . . . d 3 X 



Vn 



(V* - V) Nt oc ( 1 - 



V T 



V_ 



N T 



1 - 



v_ 



N T 



V 



N T 



exp 



PoV 



(16) 



(tn t / Po )j ~ r v r 

where it has been assumed, Vr — > oo, iVr — > oo, VyVr — >■ 0, and the ideal gas equation for 
the thermostat pressure, p = TN T /Vr, has been used. 
One obtains from Eqs. 



/v(xi, • • • ,Xjv, Pi, . . . ,p 



Ni 



1 

— exp 



p V + E 



(17) 



Z pce (p ,T,iV) N\ V T 
where the system energy E in (TTTj) depends on particle coordinates and momenta according to 
Eq. ([3]). The function Z pce (p ,T, N) in Eq. fflTj) is the (po,T,N) Pressure Canonical Ensemble 
(PCE) partition function. It is defined by the normalization condition, 



N r 



(2tt) 3 



/v(xi,...,xjv; pi, 



and it equals to, 



Zp C eG°0,^, AT) 



exp(-^) Z ce (V,T,N) 



(19) 



where Z ce (V,T, N) is the CE partition function given by Eq. (J2]). 

Using Z ce = exp(— F/T), the probability volume distribution in the PCE can be presented 

as: 



W pce (V) 



1 



exp 



F(V,T,N) + ppV 
T 



(20) 



Zpce(P0,T, N) 

First, one finds the maximum of W pce (V), which defines the most probable value of the volume, 
V = V , 

L uv Jy=v 

By definition, — dF/dV is the CE pressure flS]). Thus, the PCE equation of state fl2TT) has clear 
physical meaning: the internal pressure p for the most probable volume Vo equals to the fixed 
external pressure p , 



p(V ,N,T) = p . 



(22) 
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In the TL the volume distribution ( 1201) can be approximated as 



W pce {V) 



x exp 



1 



exp 



Zpce(P0,T, N) 

1 fd 2 F(V,T,N) 



F(V ,T,N) + Po V 



2T 



(V 



T 



Vn 



(23) 



V=V 



= (2-kujv V 



-1/2 



exp 



(y - v y 

2lu v Vq 



Thus, the most probable volume, Vo, and the average volume V are equal to each other in the 
TL, 



V 



dV V W pce {V) V 



(24) 



and uiy introduced in Eq. fl23l) defines the scaled variance of the volume fluctuations, 



V 2 - v 2 



ujv = — 

V 



T 
Vo 



r d 2 F(V,T, N) 
dV 2 



i -l 



V=V 



For the ideal gas pressure (jHJ) one gets: 

V = V 



NT 
Po 



T 
Vo 



T 
Po 



dp(V, T, N) 
dV 



(25) 



V=V 



(26) 



The Eqs. ( )23ll25l) are valid if (dp/dV) is negative. This is the case for the 'normal' equation of 
state, but is not valid for either a 1st order phase transition or critical point. Note that the TL 
N, Vq — > oo, with N/Vq being finite, is assumed to make clear notions of the phase transition 
or critical point. 



1st Order Phase Transition. In the casei of the 1st order phase transition dp/dV = 19] for 
V\ < V < V2 in the mixed phase. According to Eq. (1201) this leads to W pce (V) = const = 
(y 2 — Vi) _1 for Vi < V < Vi- Introducing the notation, V% = jVi, with 7 = const > 1, one 
finds at V\ — > 00, 



V 



V 2 



1 

V 



V W pce (V)dV 

v 2 



7 + 1 



Vl 



V 2 W pce (V)dV - V 



7 



7 + 1 



V . 



(27) 
(28) 



Thus, the volume fluctuations are anomalously large, u v oc V — >• 00. 



Critical Point. At the critical point it follows [1 

d 2 p 



dp 
W 



dV 2 



(29) 
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In this case, Eq. (1201) takes the form 



W pce {V) 



x exp 



1 



exp 



Z pce (Po,T, N) 
1 1 fd 4 F(V,T,N) 



F{V ,T,N) + Po Vo 



T 4! 



dV 4 



T 

(V - V ) 4 



(30) 



V=V 



2 A 



1/4 



r(i/4) 



exp [ - a (y - Vo) 4 ] 



where A = — (24T) 1 (9 3 p/c?V 3 )v-=v r o- The volume distribution fl30l leads to the scaled variance 



for the volume fluctuations 

v*-v 2 

OJy = 



i r(3/4) 



(31) 



where the Gamma functions in Eqs. 
one finds for (13"IT). 



are T(l/4) ^ 3.626 and T(3/4) ^ 1.225. Finally, 



1.656 T 1 / 2 



<9 3 p 

m/ 3 



y=vb 



'1/2 



(32) 



In order to estimate the scaled variance fl32l of the volume fluctuations at the critical point we 
use the van der Waals (VdW) equation of state ljj, 

N 2 \ 



p + a 



V 2 J 



(V - 6 JV) = iVT . 



(33) 



The critical point is defined by the following equations 

dp NT 2N 2 a 



dV 
d 2 p 

W 2 = (V - Nbf 



(V-Nb) 2 + V 3 
2NT 6N 2 a 



0, 



V A 



(34) 
(35) 



They give, 



T 



8 a 
27 b 



V cr = 3Nb, p c 



1 a 
27 62 



(36) 



At the critical point one finds, 

d 3 p 
dV 3 



6NT 



(v - my 



+ 



24N 2 a 



a 1 



9 2 b 5 N 3 



(37) 



Thus, for T = T cr , p = p cr , and V = V = V cr , and Eqs. (I32II37I) give, 



1.656 T 1 ! 2 



Q3p 

W 3 



1-1/2 



y=v cr . 



2.703 6 ViV . 



(3* 
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The volume fluctuations at the critical point are anomalously large, ay oc N 1 / 2 ocV — > oo. 

Similarly to the PCE we introduce now the (sq, E, N) Pressure Micro-Canonical ensemble 
(PMC), 

POO 

Z pmc (s ,E,N) = / dVex P (-s V)Z mce (V,E,N) , (39) 
Jo 

where the parameter so — Vo/ To is defined by the external conditions of the thermostat pressure 
Po and the thermostat temperature To. The micro-canonical energy E is assumed to be fixed. 
Thus, there is no energy exchange between the 'thermostat' and considered system, and the 
internal MCE temperature T (fT2j) may differ from T . Using Eq. (flOl) . the volume distribution 
in the PMC can be presented in the form, 

W pmc (V) = E N) ex P [- s V + S(V,E,N)} (40) 



1 

exp 



s Vo + S(V ,E,N) + ~ (0) (V-K>) 2 



(s ,E,N) 

The most probable (and average) volume Vq is defined by the condition, 

- s + 
The Eq. f )4T|) corresponds to, 



(42) 

which is similar to Eq. (|22|) in the PCE. The scaled variance of the volume fluctuations uoy in 
the PMC equals to: 



uj v = [V (d 2 S/dV 2 )v=v ] 1 = - tt 



T 

o 



-I -1 



(43) 

V=V 



which is again similar to Eq. (1251) in the PCE. Thus, the average volume Vo and scaled variance 
of the volume fluctuations uy are identical in the PCE and PMC in the TL. For the ideal gas 
one finds, 

N 1 

(V) pmc = Vo = - , CUy = - , (44) 
So s 

which coincide with fl26l) in the PCE. 
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C. Grand Canonical Ensemble 



The (V, T, jj) Grand Canonical Ensemble (GCE) partition function of non-relativistic non- 
interacting Boltzmann particles reads: 

Zgjy^y) = J2 exp (V) Z °*( V > T > N )> ( 45 ) 

N=0 ^ ' 

where Z ce (V,T, N) is given by Eq. (|2J) and // is the chemical potential. The GCE pressure is 
calculated as, 

p(T,fi) = I ln[Z gce (V,T,/i)] , (46) 

and it does not depend on the volume V in the TL. The CE, MCE, and GCE are thermo- 
dynamically equivalent at V — > oo. Note that the GCE is the most convenient one from the 
technical point of view. 

There is, however, an evident problem in the formulation of the (po,T,fi) Pressure Grand 
Canonical ensemble (PGC). The PGC partition function is obtained by extending Eq. ( jl~9l) . 

Z P g C (po,T,fi) = dV exp (- Z gce (V,T,fi) = — - ^— ^ , (47) 

where p(T, /i) is the GCE pressure ( l4*6|) . The (po, T, /i)-ensemble has a unique property. Among 
8 possible ensembles this is the only one where the system description includes only intensive 
quantites, Po,T and \i. For p(T, //) = p , the system volume is undefined. In the domain 
p(T,fi) > po, the PGC partition function does not exist as the integral over the volume in 
Eq. ( 1471) diverges. For p(T,p) < po, the volume distribution has the form: 

1 

Z pgc (p ,T, fj) 



W pgc (V) = - ( - - exp 



- V 



Po - P 



T 

The most probable volume Vq equals to zero, but the average volume is: 



(48) 



/* oo rp 

V = / dVV W pgc (V) = — - . (49) 

Jo Po ~ p(T,fi) 

This special ensemble will be discussed further in details for the ultra-relativistic gas in the 
next section. 
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III. ULTRA- RELATIVISTIC GAS 



In this section several examples of the pressure ensembles for the ultra-relativistic (m = 0) 
ideal gas of Boltzmann particles are considered. For simplicity only statistical systems without 
conserved charges are discussed. Thus, the number of particles is not restricted and chemical 
potential equals to zero. Furthermore, the Boltzmann statistics is used and the degeneracy 
factor is assumed to be one. 



A. Grand Canonical and Grand Micro- Canonical Ensembles 

The (V, T) GCE 4 partition function of massless non-interacting neutral Boltzmann particles 
reads: 



*~<™ = Es ^ ) L n**>~* (-1) = e m ( 

N=® v / JU i=l Af=0 v 



VT ? - 



exp(iV) ; (50) 



where N = (N) gce = VT Jir is the GCE average number of particles. The GCE system 



pressure and average energy are: 



rp rpA 

p = — lnZ gce = — = T n(T) , (51) 

V 7T Z 



(E) gce = E = T* = VT 4 = e(T) V , (52) 

where n(T) = N/V = T 3 /tt 2 and e(T) = 3T 4 /tt 2 are the particle number density and energy 
density, respectively. The GCE multiplicity distribution has the Poisson form, 

P gce (N;V,T) = — exp(-F) , (53) 
and the scaled variance of particle number distribution (153]) equals to: 

^gce = 7T77 — = 1 • (54) 



(N) 



yet 



The chemical potential connected to the number of particles equals to zero. 
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The (V, E) Grand Micro-Canonical ensemble 5 (GMC) partition function is 

00 1 / y \ iV r w / N \ 

z gmc (v,E) = 53 _ — / h^sIe - j> 

N=l ' V ' Jo i=l \ j=l / 

00 1 ™ ^4 / 4 5 4\ 

s Xy^) = lElpFTTB = 2E °M ; 3'3' 2i 27) ■ (55) 



where 0-F3 is the generalized hyper-geometric function 21[, and 



VE 3 

A = -jt ■ (56) 



7T 2 



The GMC particle number distribution function equals to 

W N (V,E) 1 

Z gmc (V,E) ~ Z gmc (V,E) EM (3JV-1)! 



P flmc (iV; V, £) = ; KAf V T ; ^ = „ 7T7 ^ 7\\ ■ (57) 



n 

It is defined for N > 1. The average number of particles in the GMC equals to [8j: 



.4 



1/4 



= ( 27 ) • (58) 

In the large volume limit the mean multiplicities in the GME and the GCE are equal, (N) gmc = 
N providing E = E and the GMC and the GCE volumes are equal. The temperature and the 
pressure in the GME are equal to: 

E fu 2 E\ 1/A ~NT T 4 

V = -FT- = -T > (59) 



37V \3V J ^ V 7T 2 

and they coincide with the corresponding quantities ( 15BI521) in the GCE. 

For N ^> 1 the particle number distribution in the GME ( 1571) can be approximated by the 
Gaussian one: 

(N - N) 



P gmc (N;V,E) S (2nco gce N) 1/2 exp 
with the GME scaled variance 



2 UJ gmc N 



21 

(60) 



(N 2 ) gmc - (N)l mc 1 



Ugmc = 777^ = 7 • (61) 



The Poisson distribution P gce (N) (1531 for N ^> 1 can be also approximated by the Gauss 
distribution (IBTil) . but with o; ffce = 1. In the TL the particle number distributions in the 
GME and GCE have both the Gauss form [l^]. The average number of particles is the same, 
(N) gmc = N, but the scaled variance is different, u gce = 1 and 0J gmc = 1/4. 



5 The energy is fixed, but the number of particles is not. The chemical potential connected to the number of 
particles equals to zero. Thus, this ensemble is named the Grand Micro-Canonical Ensemble 
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B. Pressure Grand Canonical Ensemble 



The (po, T) Pressure Grand Canonical ensemble (PGC) is defined as the ensemble with the 
fixed external pressure p$, temperature T, and the chemical potential connected to the number 
of particles equals to zero. The PGC partition function equal to [20! ]: 

Z pgc (Po,T) = J~dV exp(- ^) Z gce (V,T) = - _^ (T) , (62) 

where the relation (loTi) between p(T) and Z gce has been used. The value of p has a physical 
meaning of the external pressure. A convergence of the integral over the volume in Eq. (I6"2l 
requires the inequality, 

p(T) = \ T 4 < p . (63) 



TV 2 



Thus, at each fixed value of po, there is a 'limiting temperature' T* in the PGC: 

T < T* = (vr 2 Po ) 1/4 • (64) 
The probability volume distribution in the PGC is: 



It gives the most probable volume equal to zero, V = 0, and the average value, 



(65) 



(V) P9C = V = [°°dVV W pgc (V) = - T ^y?o,T) = T 

Jo d Po Po ~ P[T) 

which is not equal to zero and may even go to infinity at T — > T*. Using Eq. (166|) . the volume 
distribution fl65|) can be written as, 

W pgc (V) = V' 1 exp (- V/V) . (67) 

For the mean values of energy and particle number in the PGC one finds: 

{E) pgc = / dV W pgc (V) (E) gce = — T 4 V = s(T) V, (68) 
Jo 7r 

r°° 1 _ _ 

(N) pgc = / dV W pgc (V) (N) gce = — T 3 V = n(T) V , (69) 
Jo 71 

where n(T) and e(T) are the particle number density and energy density of the GCE given by 

Eq. ([ST]) and Eq. fl52l . respectively. 
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The independent variables are (V,T) in the GCE and (po,T) in the PGC If the GCE 
volume is chosen to be equal to the average volume of the PGC, then Eqs. ( 1681) and ( 1691) give: 
(E)p g c = E and (N) pgc = N, i.e. the GCE and PGC are thermodynamically equivalent. The 
fluctuations of E and N are however different in these two ensembles. Calculating, 



(V\ gc 
one finds, 



dV V 2 W pgc (V) 



T 2 d 2 Z, 



/><!'' 



Zpgc &Po 



T 



Po ~ P(T) 



2V 



(70) 



OJy 



(V 2 ) pgc 



(V)pgc 



V . 



(71) 



(V) Pgc 

Thus, in the TL limit V — > oo the volume fluctuations become anomalously large. This, in turn, 
leads to anomalous energy and particle number fluctuations. The particle number distribution 
in the PGC has the form of the geometrical distribution: 



POO 

P pgc (N; Po ,T) = / dVW pgc {V) P gce {N,V,T) = {l-r))r] N 
Jo 



(72) 



where P gce (N, V,T) has been taken from Eq. fl53|) . and rj = (T/T*) 4 < 1 . The most probable 
number of particles is iV = 0, whereas the average value (169|) is larger than zero. It equals to 
(N) pgc = — rj)^ 1 and may even go to infinity at rj — ► 1. This happens if T — > T*. From 
Eq. (G2D it follows: 



(N 2 ) pgc - (N) 



p<jc 



1 + 



V 



1 + n(T) V . 



(73) 



(N) Pg c " ' l-V 

The first term in the r.h.s. of Eq. ( 1731) corresponds to the Poisson fluctuations ( 15^1) of the GCE 
at fixed volume, whereas the second term comes from the volume fluctuations at fixed particle 
number density. The multiplicity distribution (1721) can be rewritten as, 



P pgc {N;p ,T) = P pgc (N;N) 



For iV 3> 1 the distribution P pgc approches: 



1 



N+l 



exp 



- N In 1 



1 

w 



(74) 



P pgc (N;N) 



^ 1 



exp — 



N 



(75) 



The particle number distribution P pgc (N) ([75]) satisfies the so called KNO-scaling [l^ |. 
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C. Pressure Grand Micro-Canonical Ensemble 



The (sq, E) Pressure Grand Micro-canonical ensemble (PGM) is defined as the statistical 
ensemble with fixed energy E and fixed external parameter sq = Po/T . As before, the chemical 
potential connected to the number of particles equals to zero. The PGM partition function is 
equal to: 



Zpgrn^Oi 



poo 

E) = dV exp(- s V) Z gmc (V,E) 
Jo 



(76) 



1_ ~ B N 



where 



B 



E 3 
s vr 2 



(77) 



Using Eq. [22j , 



E 



k=0 



(3k)\ 



exp(x) + 2 exp(— x/2) cos ( x 



(78) 



one finds at B > 1, 



ln[Z pgm (s ,E)} * B 1 ' 3 . 



(79) 



In the TL, the volume distribution in the PGM can be found from ( 1761) using the asymptotic 
behavior of 0-F3 function 21], 

fVE 3 \ 1/A 



\n[Z gmc (V,E)] = In 



A 
2E 



4 5 A 



0F3 I ; -, - , 2 



' 3' 3' '27 



The volume distribution in the PGM is then proportional to, 



WpgmiV) oc exp 



VE l 



3\ 1/4 



27tt 5 



soV 



exp 



0(^0 



1 / d 



2 V dV 2 



= exp[0(y)] 
(V - K,) 2 



0) 



V=Vb 



The most probable volume Vq in the PGM is defined by the condition, 

(E 3 ^ /4 
I or 

V=V 



30 

dV 



27tt 2 V 3 



- s = . 



rs2) 
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This gives, 



v ° = (83) 



One also finds, 



Thus, the volume distribution in the PGM can be approximated as: 



34) 



W pgm (V) = (2nuj v V y 1/2 



exp 



(V - Vo) 

2 UJy Vo 



55) 



where ujy = 4/ (3sq) is the scaled variance of the volume fluctuations in the PGM. The average 
volume and its fluctuations in the PGM can be also calculated in terms of the partition function 
Z pgm (s ,E) using Eqs. ((ZEHZ2J), 



(V) pgm = I dVW pgm (V) V = - s j-A^ , (86) 



/■oo i o2 7 7712 

(O— = / "VW—W ^ = J- - —^875 . (87) 



pgm r^_, 



(V/pgm 3So 

The condition (|83|) can be written as, 

*o = ~ 2 T 3 = I , (89) 

and it means that the internal pressure p equals to sqT. The most probable volume Vo ( 1831) and 
average volume (^)psm fl86|) are then equal to each other, and both are equal to En 2 / (3T 4 ). This 
corresponds to the fixed volume in the MCE with fixed energy E and the MCE temperature 
T. The volume distribution in the PGM is therefore different from that in the PGC (l67j) . In 
contrast to the PGC, there is an extensive variable, the energy E, in the PMG. This leads to the 
system average volume proportional to the energy and given by Eq. (|83p . The internal pressure 
equals to p — SqT fl89|) . As a result, the volume fluctuations in the PGM are Gaussian (1831) 
with the finite scaled variance (188]) . Using Eq. (J89]) the scaled variance (l8"8"j) of the volume 
fluctuations can be expressed in terms of the MCE temperature, ioy = An 2 /(3T 3 ). One also 
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finds: 



<iV> 



pgm 



(N 2 ) 



pgm 



'-'pgm 
1 

Zpgm 

(N 2 ) 



TD dZ pgm 1 j/g 

B ^B~ = 3 B > 

d OZpgm ^ 1 Wg 

~dB~ = 9 * 



9 



^2/3 



pgm 



- w 



pgm ^ 



(iV)p 3 m " 3 • 

In the TL, the multiplicity distribution in the PGM can be approximated as: 



P pgm (N- s , E) S (2 



7T 



N) 



-1/2 



cxp 
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(N-N) 

2 ^Jpgm -A" 



(90) 
(91) 
(92) 

(93) 



J B 1/3 /3 = (N) pgrn and cu p9m 



1/3 . The scaled variance in the PGM can be 



(94) 



where iV 

presented in the TL as the following, 

1 111 

Upgm — ^mce + n ~ ^ + ^2 _ 3' 

where n = T 3 /n 2 is the MCE particle number density. The first term in the r.h.s. of Eq. (|9l 
Umce = 1/4 , is due to the particle number fluctuations in the MCE with fixed volume, and the 
second term is the contribution due to the volume fluctuations. 



IV. ENSEMBLES WITH EXTERNAL VOLUME FLUCTUATIONS 



The multiplicity distributions P(N) in relativistic gases [7|, |8|, |9j, [10J, lUl, [12|, [13] are sensitive 
to conservation laws obeyed by the system, and therefore to fluctuations of extensive quantities 
E and Q. The examples considered in the previous section demonstrate that the volume fluctu- 
ations also influence the particle number fluctuations. Thus, for the calculation of multiplicity 
distributions, the choice of the statistical ensemble is then not a matter of convenience, but a 
physical question. On the other hand, the fluctuations of extensive quantities A = (V, E, Q) 
depend not on the system's physical properties, but rather on external conditions. One can 
imagine a huge variety of these conditions, thus, 8 statistical ensembles discussed in the previous 
sections are only some special examples. 

A more general concept of the statistical ensembles based on Eq. (CQ) was recently suggested 



in Ref. 



161 ] . The system volume may exhibit fluctuations described by the externally given 



distribution. When V is the only fluctuating variable, Eq. ([T]) is reduced to 

P a (0) = [dV P a (E, V, Q) P mce (0; E, V, Q) , 



(95) 
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P a (V) = (2ttu v a 2 v V) 1/2 exp 



.21 



(96) 



where P a {V) is externally given volume distribution. 

The effect of volume fluctuations on the particle number fluctuations is calculated for 
the system of non-interacting massless Boltzmann particles with zero chemical potential. At 
fixed volume the system is treated within the GMC, and the particle number distribution is 
P gmc (N;V,E) ([57D. 

In the first example, the volume distribution is assumed to be: 

2(jJ V CtyV 

where Uy = 4n 2 / (3T 3 ) coincides with (1551) in the PGM. The choice of P a (V) results in a simple 
correspondence to the GMC and PGM in the TL. In Eq. (1961) . ay is a dimensionless tuneable 
parameter which determines the width of the distribution. In the limit ay — > 0, Eq. (1961) 
becomes a Dirac 5-function, 5(V — V). This corresponds to the GMC. For ay = 1, Eq. fl96l) 
results in the PGM volume fluctuations ([85]) in the TL. The particle number distribution reads, 

POO 

P a (N) = / dV P a {V) P gmc {N;V,E) . (97) 
Jo 

A substitution of P a (V) in Eq. (J9"7j) by the distribution fl96|) results in P a {N) in the TL given 
by the Gaussian, 

(N - N)' 



P a (N) = (2iru a N) 1/2 exp 



2u n N 



(98) 



where the average number of particles iV is defined by the energy and average volume, 
N=[V E 3 /(277r 2 )] 1 / 4 , and the scaled variance of the particle number distribution fl98|) equals, 

fa = ^ + ^ a 2 v . (99) 

The first term in the r.h.s. of Eq. (199|) equals to the GMC scaled variance at fixed E and V, the 
second term is due to the volume fluctuations. As it can be expected, u a = u gmc for ay = 0, 
and uj a = 0J pgm for ay = 1. It also follows, uj gmc < u a < u> pgm for < a v < 1, and u a > ui pgm for 
ay > 1. The a-ensemble defined by Eqs. f l97] [961) presents an extension of the GMC (ay = 0) 
and PGM (ay = 1) to a more general volume distribution. 



In the second example both E and V are assumed to fluctuate. Equation fl97|) should be 
then extended as: 

P a (N) = J dEdV P a (E, V) P gmc (N; E, V) . (100) 
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First, uncorrelated volume and energy distributions are considered: 



P a (E,V) = P 1 (V)xP 2 (E), (101) 
where -Pi(V) is given by Eq. ( 1961) and P2{E) is taken in the following form, 



P 2 (E) = (2ttu e a\E) 1/2 exp 



(E - E)' 

2uj e a% E 



(102) 



where E = 3T A V/n 2 and uo E = 4T. A substitution of P a (E, V) in Eq. (11001) by the above distri- 
butions results in P a (N) in the TL given again by the Gaussian ( |98l) with the average number 
of particles N defined by the average energy and average volume, N = [V E /(277T 2 )] 1 / 4 , and 



the scaled variance equal to: 



u a = J + ^ o| + ^Oy . (103) 



The first term in the r.h.s. of Eq. (11031) is the MCE scaled variance at fixed E and V, the second 
term is due to energy fluctuations, and the third one comes from volume fluctuations. The 
results of the GMC, PGM, and GCE for the scaled variance of the particle number distribution 
can be obtained from Eq. (11031) at a v = a E = 0; a v — 1, a E — 0; and ay = 0, a E = 1, 
respectively 

Examples of correlated V and E distributions can be constructed as follows. One assumes 
the factorized form (1 1011) with Eq. (11021) for the energy distribution, but with E = 3T a V/tt 2 
and uj e = 4T, i.e. the average energy E depends on the volume V. Assuming, 

Pi(V) = i= exp(- V/V) , (104) 

with V given by Eq. (|66j) . from Eqs. (jlOip and (j!02p with a E = 1, the results for the PGC are 
obtained. 

The relation, E = 3T a V/tt 2 , can be used together with any form of the volume distribution 
P\{V). This yields a generalization of the GCE to the systems with externally given volume 
fluctuations. The energy fluctuations at fixed V can be also selected in accordance with the 
physics requirements. 

As the third example of the statistical ensembles with volume fluctuations one refers to the 
recent paper 17j where the micro-canonical ensemble with the scaling volume fluctuations for 



the ultra-relativistic ideal gas has been considered. The volume fluctuations were assumed to 
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have the specific scaling properties. They were chosen to describe the KNO scaling [18| of 
the particle multiplicity distributions measured in proton-proton collisions at high energies. 
A striking feature of the model is power law form of the single momentum spectrum at high 
momenta, instead of the exponential Boltzmann distribution in the systems with fixed volume. 

V. SUMMARY 

In this paper the volume fluctuations in the statistical mechanics have been studied. The 
statistical ensembles with fixed external pressure were considered. Statistical systems of classi- 
cal non-relativistic particles and non-interacting massless particles were discussed. The volume 
fluctuations in the pressure canonical ensemble become anomalously large in the thermody- 
namic limit for the first order phase transition, ojy oc N, or at the critical point, lo oc A 1 / 2 . 
Another type of anomalous volume fluctuations takes place for the pressure grand canonical 
ensemble. In this special ensemble, all thermodynamical variables are the intensive quanti- 
ties, (po,T,fi). In the thermodynamical limit the mean particle multiplicity obtained within 
the considered ensembles with volume fluctuations is equal to the mean calculated within the 
micro-canonical, canonical, and grand canonical ensembles. This is not valid, however, for the 
scaled variance. The influence of the volume fluctuations in the pressure ensembles on particle 
number fluctuations have been discussed for ultra-relativistic ideal gas. 

Following Ref. [16J the ensembles with volume fluctuations determined by the externally 
given distribution function were introduced. Multiplicity fluctuations are sensitive to the vol- 

he volume fluctuations may also influence the behavior of the particle 



ume fluctuations. r . 

momentum spectra 12J. Thus, we believe that the concept of statistical ensembles with fluc- 
tuating extensive quantities, in particular, with fluctuating volume of the statistical system, 
may be appropriate for the statistical description of hadron production in relativistic collisions. 
It may be also useful for other physical systems. In fact, in all cases when the equilibrium 
statistical mechanics is used to calculate the fluctuations of the system properties. 
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